Protein-coding Genes Research Articles

Plastome structure, phylogeny and evolution of plastid genes in Reevesia (Helicteroideae, Malvaceae).

Reevesia is an eastern Asian-eastern North American disjunction genus in the family Malvaceae s.l. and comprises approximately 25 species. The relationships within the genus are not well understood. Here, 15 plastomes representing 12 Reevesia species were compared, with the aim of better understanding the species circumscription and phylogenetic relationships within the genus and among genera in the family Malvaceae s.l. The 11 newly sequenced plastomes range between 161,532 and 161, 945bp in length. The genomes contain 114 unique genes, 18 of which are duplicated in the inverted repeats (IRs). Gene content of these plastomes is nearly identical. All the protein-coding genes are under purifying selection in the Reevesia plastomes compared. The top ten hypervariable regions, SSRs, and the long repeats identified are potential molecular markers for future population genetic and phylogenetic studies. Phylogenetic analysis based on the whole plastomes confirmed the monophyly of Reevesia and a close relationship with Durio (traditional Bombacaceae) in subfamily Helicteroideae, but not with the morphologically similar genera Pterospermum and Sterculia (both of traditional Sterculiaceae). Phylogenetic relationships within Reevesia suggested that two species, R. pubescens and R. thyrsoidea, as newly defined, are not monophyletic. Six taxa, R. membranacea, R. xuefengensis, R. botingensis, R. lofouensis, R. longipetiolata and R. pycnantha, are suggested to be recognized.

LncRNAway: a web-based sgRNA design tool for precise and effective suppression of long noncoding RNAs.

Thousands of long noncoding RNAs (lncRNAs) have been annotated via high-throughput RNA sequencing, yet only a small fraction have been functionally investigated. Genomic knockout is the mainstream strategy for studying the biological function of protein-coding genes and lncRNAs, whereas the complexity of the lncRNA locus, especially the natural antisense lncRNAs (NAT-lncRNAs), presents great challenges. Knocking out lncRNAs often results in unintended disruptions of neighboring protein-coding genes and small RNAs, leading to ambiguity in observing phenotypes and interpreting biological function. To address this issue, we launched LncRNAway, a user-friendly web tool based on the BESST (branchpoint to 3' splicing site targeting) method, to design sgRNAs for lncRNA knockout. LncRNAway not only provides specific and effective lncRNA knockout guidelines but also integrates genotyping primers and quantitative PCR primers designing, thereby streamlining experimental procedures of lncRNA function study. LncRNAway is freely available at https://www.lncrnaway.com.

Intraspecific and interspecific variations in the synonymous codon usage in mitochondrial genomes of 8 pleurotus strains

In this study, we investigated the codon bias of twelve mitochondrial core protein coding genes (PCGs) in eight Pleurotus strains, two of which are from the same species. The results revealed that the codons of all Pleurotus strains had a preference for ending in A/T. Furthermore, the correlation between codon base compositions and codon adaptation index (CAI), codon bias index (CBI) and frequency of optimal codons (FOP) indices was also detected, implying the influence of base composition on codon bias. The two P. ostreatus species were found to have differences in various base bias indicators. The average effective number of codons (ENC) of mitochondrial core PCGs of Pleurotus was found to be less than 35, indicating strong codon preference of mitochondrial core PCGs of Pleurotus. The neutrality plot analysis and PR2-Bias plot analysis further suggested that natural selection plays an important role in Pleurotus codon bias. Additionally, six to ten optimal codons (ΔRSCU > 0.08 and RSCU > 1) were identified in eight Pleurotus strains, with UGU and ACU being the most widely used optimal codons in Pleurotus. Finally, based on the combined mitochondrial sequence and RSCU value, the genetic relationship between different Pleurotus strains was deduced, showing large variations between them. This research has improved our understanding of synonymous codon usage characteristics and evolution of this important fungal group.

The first high-quality genome assembly and annotation of Lantana camara, an important ornamental plant and a major invasive species

This study presents the first annotated, haplotype-resolved, chromosome-scale genome of Lantana camara, a flowering shrub native to Central America and known for its dual role as an ornamental plant and an invasive species. Despite its widespread cultivation and ecological impact, the lack of a high-quality genome has hindered the investigation of traits of both ornamental and invasive. This research bridges the gap in genomic resources for L. camara, which is crucial for both ornamental breeding programs and invasive species management. Whole-genome and transcriptome sequencing were utilized to elucidate the genetic complexity of a diploid L. camara breeding line UF-T48. The genome was assembled de novo using HiFi and Hi-C reads, resulting in two phased genome assemblies with high Benchmarking Universal Single-Copy Orthologs (BUSCO) scores of 97.7%, indicating their quality. All 22 chromosomes were assembled with pseudochromosomes averaging 117 Mb. The assemblies revealed 29 telomeres and an extensive presence of repetitive sequences, primarily long terminal repeat transposable elements. The genome annotation identified 83,775 protein-coding genes, with 83% functionally annotated. In particular, the study mapped 42 anthocyanin and carotenoid candidate gene clusters and 12 herbicide target genes to the assembly, identifying 38 genes spread across the genome that are integral to flower color development and 53 genes for herbicide targeting in L. camara. This comprehensive genomic study not only enhances the understanding of L. camara’s genetic makeup but also sets a precedent for genomic research in the Verbenaceae family, offering a foundation for future studies in plant genetics, conservation, and breeding.

A chromosome-level genome assembly of an avivorous bat species (Nyctalus aviator)

Currently, three carnivorous bat species, namely Ia io, Nyctalus lasiopterus, and Nyctalus aviator, are known to actively prey on seasonal migratory birds (hereinafter referred to as “avivorous bats”). However, the absence of reference genomes impedes a thorough comprehension of the molecular adaptations of avivorous bat species. Herein, we present the high-quality chromosome-scale reference genome of N. aviator based on PacBio subreads, DNBSEQ short-reads and Hi-C sequencing data. The genome assembly size of N. aviator is 1.77 Gb, with a scaffold N50 of 102 Mb, of which 99.8% assembly was anchored into 21 pseudo-chromosomes. After masking 635.1 Mb repetitive sequences, a total of 19,412 protein-coding genes were identified, of which 99.3% were functionally annotated. The genome assembly and gene prediction reached 96.1% and 96.1% completeness of Benchmarking Universal Single-Copy Orthologs (BUSCO), respectively. This chromosome-level reference genome of N. aviator fills a gap in the existing information on the genomes of carnivorous bats, especially avivorous ones, and will be valuable for mechanism of adaptations to dietary niche expansion in bat species.

Complete Plastid Genome Sequences of Three Tropical African Beilschmiediineae Trees (Lauraceae: Crytocaryeae)

Millions of years of isolation have given Madagascar a unique flora that still reflects some of its relationship with the continents of Africa and India. Here, the complete chloroplast sequence of Beilschmiedia moratii, a tropical tree in Madagascar, was determined. The plastome, with a length of 158,410 bp, was 143 bp and 187 bp smaller than those of two closely related species, B. pierreana and Potameia microphylla, in sub-Saharan Africa and Madagascar with published sequences, respectively. A total of 124 repeats and 114 simple sequence repeats (SSRs) were detected in the plastome of B. moratii. Six highly variable regions, including ndhF, ndhF-rpl32, trnC-petN, pebE-petL, rpl32-trnL, and ycf1, among the three African species were identified and 1151 mutation events, including 14 SVs, 351 indels, and 786 substitutions, were accurately located. There were 634 mutation events between B. moratii and P. microphylla with a mean nucleotide variability (π) value of 0.00279, while there were 827 mutation events between B. moratii and B. pierreana with a mean π value of 0.00385. The Ka/Ks ratios of 86 protein-coding genes in the three African species were less than 1, and the mean value between B. moratii and P. microphylla was 0.184, while the mean value between B. moratii and B. pierreana was 0.286. In this study, the plastid genomes of the three African Beilschmiediineae species were compared for the first time and revealed that B. moratii and P. microphylla from Madagascar were relatively conserved, with low mutation rates and slower evolutionary rates.

The adaptive prowess of the lesser mealworm Alphitobius diaperinus (Coleoptera: Tenebrionidae): Insights from mouthpart morphology and mitogenomic analysis

The lesser mealworm, Alphitobius diaperinus (Panzer, 1797), is a pervasive pest known to infest stored products and poultry facilities worldwide, posing significant challenges for pest management. Understanding its adaptability to diverse environments is crucial for devising effective control strategies. This study investigated adult mouthparts and analyzed the mitogenome of A. diaperinus using specimens obtained from a cat litter box. The examination of mouthparts revealed notable intraspecific variations, suggesting potential adaptations to different ecological conditions. Analysis of the mitogenomes unveiled nucleotide variations in protein-coding genes (PCGs), ribosomal RNA genes (rRNAs), and the control region, with all mutations in PCGs being nonsynonymous, primarily occurring at the second codon positions. Evaluation of Ka/Ks ratios indicated purifying selection acting on the PCGs of A. diaperinus. Phylogenetic analyses consistently placed the newly sequenced sample within the same clade as previously sequenced A. diaperinus samples. The presence of A. diaperinus in cat litter highlights its remarkable adaptability to diverse habitats. Furthermore, intraspecific variation in mouthpart morphology and the mitogenome suggests ongoing adaptations to various food sources and environmental niches. These findings underscore the dynamic nature of A. diaperinus adaptations and emphasize the importance of continued research for effective pest management strategies.

The genome sequence of the Brown Ash Ermine moth, Zelleria hepariella Stainton, 1849

We present a genome assembly from a male Zelleria hepariella (the Brown Ash Ermine; Arthropoda; Insecta; Lepidoptera; Yponomeutidae). The genome sequence is 428.8 megabases in span. Most of the assembly is scaffolded into 19 chromosomal pseudomolecules, including the Z sex chromosome. The mitochondrial genome has also been assembled and is 16.31 kilobases in length. Gene annotation of this assembly on Ensembl identified 15,718 protein coding genes.

Chromosomal-level genome assembly and annotation of the tropical sea cucumber Holothuria scabra

Holothuria scabra, a commercially valuable yet ecologically vulnerable tropical holothuroid, has experienced a severe decline in its wild populations, especially in China. Genomic resources are crucial for the development of effective genomic breeding projects and stock conservation strategies to restore these natural populations. Until now, a high-quality, chromosome-level reference genome for H. scabra has not been available. Here, we employed Oxford Nanopore and Hi-C sequencing technologies to assemble and annotate a high-quality, chromosome-level reference genome of H. scabra. The final genome comprised 31 scaffolds with a total length of 1.19 Gb and a scaffold N50 length of 53.52 Mb. Remarkably, 1,191.67 Mb (99.95%) of the sequences were anchored to 23 pseudo-chromosomes, with the longest one spanning 79.75 Mb. A total of 34,418 protein-coding genes were annotated in the final genome, with BUSCO analysis revealing 98.01% coverage of metazoa_odb10 genes, marking a significant improvement compared to the previous report. These chromosome-level sequences and annotations will provide an essential genomic basis for further investigation into molecular breeding and conservation management of H. scabra.

Characterization of the Complete Mitochondrial Genome of Salassa thespis (Lepidoptera: Saturniidae) and Comparison with Other Lepidoptera Species

There are many species of Lepidoptera, but few complete mitochondrial genomes of Lepidoptera have been included in databases. Here, the complete mitochondrial genome sequence of Salassa thespis was isolated and characterized. It was 15,302 bp in length and contained 13 protein-coding genes (PCGs), two rRNA genes, 22 tRNA genes and an A + T-rich region. Among the 13 PCGs, the initiation codon of cytochrome c oxidase subunit 1 (cox1) was CGA, and the rest were ATN. The cox1 and cox2 genes had an incomplete stop codon T, while the rest terminated with TAA. Codon usage analysis showed that Phe, Ile, Leu and Asn were the most frequent amino acids, while Trp was the least. Like other Lepidopterans, some conserved motifs were found in the A + T-rich region, including a 17 bp poly-T guided by ATAGA, the AT-rich area and a poly-A element. Bayesian inference and maximum likelihood phylogenetic tree analysis based on 13 PCGs of S. thespis confirmed that it belonged to the Saturniidae family and showed the following relationship: (S. thespis + (Rhodinia fugax + Samia canningi)). The enrichment of mitochondrial DNA provides reference information for the study of the evolution and diversity of Lepidoptera insects.

Mitochondrial Genome Resource of the Cottony Ash Psyllid, a Host of a Newly Identified ‘Candidatus Liberibacter’ bacterium

Cottony ash psyllid (CAP, Phyllopsis discrepans) is an important, invasive insect pest of ash trees in North America where it has established populations and is the host of a newly identified strain of 'Candidatus Liberibacter solanacearum'. However, not much is known about the diversity of its introduced population. In this study, a CAP mitochondrial genome (mitogenome) sequence was obtained from a collection in Saskatoon, Saskatchewan, Canada. The CAP mitogenome is a circular DNA of 18,824 bp, encoding 13 protein-coding genes, 21 transfer RNA genes, and 2 ribosomal RNA genes. BLAST search using the CAP mitogenome as query against GenBank sequence database showed the mitogenome of Euphyllura phillyreae (15,202 bp) was the most similar (Query Coverage = 77%; Percentage Identity = 78.90%). The CAP mitogenome is significantly different from other known psyllid mitogenomes with the presence of a 4,357 bp control region (CR). The mitogenome sequence will further genomic understand of CAP.

Transposable elements alter gene expression and may impact response to cisplatin therapy in ovarian cancer.

Cisplatin is widely employed for cancer treatment; therefore, understanding resistance to this drug is critical for therapeutic practice. While studies have delved into differential gene expression in the context of cisplatin resistance, findings remain somewhat scant. We performed a comprehensive investigation of Transposable Elements (TEs) expression and their impact in host genes in two cisplatin-treated ovarian cancer cell lines. RNA-seq, ATAC-seq, and in-depth bioinformatics analysis were used to compare cisplatin-sensitive and -resistant ovarian cancer cell lines. Our results reveal that cisplatin therapy alters the expression of protein-coding genes, but also key TEs, including LINE1, Alu, and endogenous retroviruses, in both cisplatin-sensitive and -resistant cell lines. By co-expressing with downstream genes or by creating chimeric transcripts with host genes at their insertion sites, these TEs seem to control the expression of protein-coding genes, including tumor-related genes. Our model uncovers TEs influencing the expression of cancer genes and cancer pathways. Collectively, our findings indicate that TEs alterations associated with cisplatin treatment occur in critical cancer genes and cellular pathways synergically. This research highlights the importance of considering the entire spectrum of transcribed elements in the genome, especially TE expression, for a complete understanding of complex models like cancer response to treatment.

Deciphering male influence in gynogenetic Pengze crucian carp (Carassius auratus var. pengsenensis): insights from Nanopore sequencing of structural variations

In this study, we investigate gynogenetic reproduction in Pengze Crucian Carp (Carassius auratus var. pengsenensis) using third-generation Nanopore sequencing to uncover structural variations (SVs) in offspring. Our objective was to understand the role of male genetic material in gynogenesis by examining the genomes of both parents and their offspring. We discovered a notable number of male-specific structural variations (MSSVs): 1,195 to 1,709 MSSVs in homologous offspring, accounting for approximately 0.52%–0.60% of their detected SVs, and 236 to 350 MSSVs in heterologous offspring, making up about 0.10%–0.13%. These results highlight the significant influence of male genetic material on the genetic composition of offspring, particularly in homologous pairs, challenging the traditional view of asexual reproduction. The gene annotation of MSSVs revealed their presence in critical gene regions, indicating potential functional impacts. Specifically, we found 5 MSSVs in the exonic regions of protein-coding genes in homologous offspring, suggesting possible direct effects on protein structure and function. Validation of an MSSV in the exonic region of the polyunsaturated fatty acid 5-lipoxygenase gene confirmed male genetic material transmission in some offspring. This study underscores the importance of further research on the genetic diversity and gynogenesis mechanisms, providing valuable insights for reproductive biology, aquaculture, and fostering innovation in biological research and aquaculture practices.

The genome sequence of the Winter Shade, Tortricodes alternella [Denis & Schiffermüller], 1775

We present a genome assembly from an individual male Tortricodes alternella (the Winter Shade; Arthropoda; None; Lepidoptera; Tortricidae). The genome sequence is 441.2 megabases in span. Most of the assembly is scaffolded into 22 chromosomal pseudomolecules, including the Z sex chromosome. The mitochondrial genome has also been assembled and is 15.73 kilobases in length. Gene annotation of this assembly on Ensembl identified 17,280 protein coding genes.

Macroevolutionary dynamics in micro-organisms: generalists give rise to specialists across biomes in the ubiquitous bacterial phylum Myxococcota.

Prokaryotes dominate the Tree of Life, but our understanding of the macroevolutionary processes generating this diversity is still limited. Habitat transitions are thought to be a key driver of prokaryote diversity. However, relatively little is known about how prokaryotes successfully transition and persist across environments, and how these processes might vary between biomes and lineages. Here, we investigate biome transitions and specialisation in natural populations of a focal bacterial phylum, the Myxococcota, sampled across a range of replicated soils and freshwater and marine sediments in Cornwall (UK). By targeted deep sequencing of the protein-coding gene rpoB, we found >2000 unique Myxococcota lineages, with the majority (77%) being biome specialists and with only <5% of lineages distributed across the salt barrier. Discrete character evolution models revealed that specialists in one biome rarely transitioned into specialists in another biome. Instead, evolved generalism mediated transitions between biome specialists. State-dependent diversification models found variation in speciation rates across the tree, but this variation was independent of biome association or specialisation. These results were robust to phylogenetic uncertainty, assumed amounts of unsampled diversity resulting in an incomplete phylogeny, and different levels of species delineation. Overall, our results are consistent with a "jack-of-all-trades" trade-off where generalists suffer a cost in any individual environment, resulting in rapid evolution of niche specialists and shed light on how bacteria could transition between biomes.

Senataxin deficiency disrupts proteostasis through nucleolar ncRNA-driven protein aggregation.

Senataxin is an evolutionarily conserved RNA-DNA helicase involved in DNA repair and transcription termination that is associated with human neurodegenerative disorders. Here, we investigated whether Senataxin loss affects protein homeostasis based on previous work showing R-loop-driven accumulation of DNA damage and protein aggregates in human cells. We find that Senataxin loss results in the accumulation of insoluble proteins, including many factors known to be prone to aggregation in neurodegenerative disorders. These aggregates are located primarily in the nucleolus and are promoted by upregulation of non-coding RNAs expressed from the intergenic spacer region of ribosomal DNA. We also map sites of R-loop accumulation in human cells lacking Senataxin and find higher RNA-DNA hybrids within the ribosomal DNA, peri-centromeric regions, and other intergenic sites but not at annotated protein-coding genes. These findings indicate that Senataxin loss affects the solubility of the proteome through the regulation of transcription-dependent lesions in the nucleus and the nucleolus.

A high-quality chromosome-scale genome assembly of blood orange, an important pigmented sweet orange variety

Blood orange (BO) is a rare red-fleshed sweet orange (SWO) with a high anthocyanin content and is associated with numerous health-related benefits. Here, we reported a high-quality chromosome-scale genome assembly for Neixiu (NX) BO, reaching 336.63 Mb in length with contig and scaffold N50 values of 30.6 Mb. Furthermore, 96% of the assembled sequences were successfully anchored to 9 pseudo-chromosomes. The genome assembly also revealed the presence of 37.87% transposon elements and 7.64% tandem repeats, and the annotation of 30,395 protein-coding genes. A high level of genome synteny was observed between BO and SWO, further supporting their genetic similarity. The speciation event that gave rise to the Citrus species predated the duplication event found within them. The genome-wide variation between NX and SWO was also compared. This first high-quality BO genome will serve as a fundamental basis for future studies on functional genomics and genome evolution.

Transcriptomic Profiles Associated with Experimental Placebo Effects in Chronic Pain.

Gene expression networks associated with placebo effects are understudied; in this study, we identified transcriptomic profiles associated with placebo responsivity. Participants suffering from chronic pain underwent a verbal suggestion and conditioning paradigm with individually tailored thermal painful stimulations to elicit conditioned placebo effects. Participants reported pain intensity on a visual analog scale (VAS) anchored from zero = no pain to 100 = maximum imaginable pain. RNA was extracted from venous blood and RNA sequencing and validation tests were performed to identify differentially expressed genes (DEGs) associated with placebo effects, controlling for sex and level of pain. Unbiased enrichment analyses were performed to identify biological processes associated with placebo effects. Of the 10,700 protein-coding genes that passed quality control filters, 667 were found to be associated with placebo effects (FDR <0.05). Most genes (97%) upregulated were associated with larger placebo effects. The 17 top transcriptome-wide significant genes were further validated via RT-qPCR in an independent cohort of chronic pain participants. Six of them (CCDC85B, FBXL15, HAGH, PI3, SELENOM, and TNFRSF4) showed positive and significant (P < 0.05) correlation with placebo effects in the cohort. The overall DEGs were highly enriched in regulation of expression of SLITs and ROBOs (R-HSA-9010553, FDR = 1.26e-33), metabolism of RNA (R-HSA-8953854, FDR = 1.34e-30), Huntington's disease (hsa05016, FDR = 9.84e-31), and ribosome biogenesis (GO:0042254, FDR = 2.67e-15); alternations in these pathways might jeopardize the proneness to elicit placebo effects. Future studies are needed to replicate this finding and better understand the unique molecular dynamics of people who are more or less affected by pain and placebo.

A chromosome-level haplotype-resolved genome assembly of oriental tobacco budworm (Helicoverpa assulta)

Oriental tobacco budworm (Helicoverpa assulta) and cotton bollworm (Helicoverpa armigera) are two closely related species within the genus Helicoverpa. They have similar appearances and consistent damage patterns, often leading to confusion. However, the cotton bollworm is a typical polyphagous insect, while the oriental tobacco budworm belongs to the oligophagous insects. In this study, we used Nanopore, PacBio, and Illumina platforms to sequence the genome of H. assulta and used Hifiasm to create a haplotype-resolved draft genome. The Hi-C technique helped anchor 33 primary contigs to 32 chromosomes, including two sex chromosomes, Z and W. The final primary haploid genome assembly was approximately 415.19 Mb in length. BUSCO analysis revealed a high degree of completeness, with 99.0% gene coverage in this genome assembly. The repeat sequences constituted 38.39% of the genome assembly, and we annotated 17093 protein-coding genes. The high-quality genome assembly of the oriental tobacco budworm serves as a valuable genetic resource that enhances our comprehension of how they select hosts in a complex odour environment. It will also aid in developing an effective control policy.

Genomic analysis of Clostridium perfringens type D isolates from goat farms

C. perfringens type D strains are the leading cause of enterotoxaemia in ruminants such as goats, sheep, and cattle. However, there has been no prior research on the genomic characteristics of C. perfringens type D strains from various regions in China. Here, we investigated the antibiotic resistance, genomic characteristics, and phylogenetic relationship of C. perfringens type D isolates recovered from goat farms in Shaanxi, Gansu, and Ningxia provinces. The antibiotic resistance test indicated that the isolates displayed high minimum inhibitory concentration (MIC) values to sulfafurazole, whereas the other antibiotics tested, such as penicillin, enrofloxacin, and florfenicol, worked well on them. Additionally, only tetracycline resistance genes [tetA(P) and tetB(P)] were identified from the isolates. A collective of 13 toxin genes, including etx and cpe were detected among the isolates. Sequence comparison revealed that the etx and cpe genes shared high sequence identities, and they could coexist on a pCW3-like plasmid, representing a potential risk to both animal breeding and public health. Phylogenetic analysis using core genome multi-locus sequence typing (cgMLST) and core genome single nucleotide polymorphisms (SNPs) revealed the close genetic relationship and potential regional/transregional transmission of the C. perfringens type D isolates in Shaanxi and Gansu provinces. Furthermore, pan-genomic analysis suggested the functional differences at the protein-coding gene level, although isolates from the same source shared a close genetic relationship. In conclusion, this study indicated the antibiotic resistance, virulence markers, potential transregional transmission, and genomic diversity of C. perfringens type D strains from various regions in China, which could provide references for the prevention of C. perfringens foodborne diseases and further research.